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Оценивание параметра мультипликативного 
шума на изображении с помощью 
мультифрактальных показателей 


Разработан метод оценивания параметра мультипликативного шума по мультифрактальному пока- 
зателю. Исследованы характеристики разработанного метода на тестовых изображениях. 


Введение 


Проблемой большинства систем интерпретации радиолокационных или медицин- 
ских изображений является большой объем информации, которая подлежит обработке. 
Вследствие этого применяются интеллектуальные системы обработки изображений, 
содержащие несколько процедур, каждая из которых имеет свою цель. Интерпретация 
результатов обработки изображений выполняется процедурами обработки высокого 
уровня и зависит от результатов работы процедур низкого уровня. При регистрации, 
кодировании и передаче изображений возникают неизбежные помехи, которые влияют 
на результат интерпретации последних. Для улучшения качества изображения перед 
процедурами обработки высокого уровня (сегментацией, идентификацией, распозна- 
ванием) применяются методы предварительной обработки изображений. 

Большинство методов предварительной обработки изображений [1], [2] и неко- 
торые методы контурной сегментации [3], [4] предполагают, что характер шума и его 
статистические параметры известны. Во многих практических задачах нет априорных 
сведений о характере шума и его статистических параметрах либо такая информация 
труднодоступна [5]. Однако эта информация определяет качество фильтрации или 
контурной сегментации изображений. Например, завышенная оценка уровня аддитив- 
ной гауссовской помехи ведет к размыванию контуров объектов, что также отражается 
на ошибке распознавания. В некоторых случаях от системы обработки изображений 
может потребоваться работа в «слепом» режиме. Тогда параметры шума идентифи- 
цируются по наблюдаемому изображению (рис. 1). 


Ввод исходного Идентификация типа Оценка параметров шума 
зашумленного изображения шума 


Вывод обработанного Предварительная Выбор структуры и параметров 
изображения обработка метода предварительной обработки 


Рисунок 1 — Структурная схема предварительной обработки изображений 


Диапазон шумов, встречающихся при обработке изображений, достаточно ши- 
рок [6], [7]. 
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В системах, использующих естественное освещение, часто возникают простран- 
ственные изменения освещенности, которые являются причиной мультипликативной 
помехи на изображении. Она моделируется как 


Их, у) = о (у) (у), (п) 
где /(х, у)- полутоновое изображение, /,(х,у) — отражающая способность объекта, 
А(х, у) — освещенность объекта, х, у — пространственные координаты. Для того чтобы 


структура и оценки параметров системы автоматизированной идентификации шума 
были обоснованными, необходимо выбрать модель помеховой ситуации. 

Статистическая модель помеховой ситуации для мультипликативного шума имеет 
вид (1). Для подавления мультипликативного шума существуют особые методы предва- 
рительной обработки изображений, однако помеховую ситуацию оценивает специалист 
в области распознавания изображений в неавтоматизированном режиме. Очевидно, что 
неавтоматизированная идентификация шума требует от оператора квалификации 
эксперта. Неточная оценка помеховой ситуации или параметров шума могут привести 
к некачественной предварительной обработке. Это обусловлено тем, что при выборе 
метода предварительной обработки и его параметров необходима априорная информация 
о процессе формирования изображений. Часто такая информация отсутствует, поэтому 
параметры метода предварительной обработки изображений подбираются путем на- 
стройки, которая не всегда позволяет отразить помеховую ситуацию. Это влияет на 
дальнейший процесс обработки изображений. Например, недооценка уровня мульти- 
пликативной помехи вызывает большие разрывы контуров объектов изображений 
при решении задачи контурной сегментации. Это, в свою очередь, влияет на ошибку 
распознавания этих объектов. Поэтому актуальной является проблема автоматизации 
процедуры идентификации шума на изображениях. 

Целью данной работы являлось снижение погрешности оценивания отношения 
сигнал/шум для мультипликативного шума на изображениях путем использования 
мультифрактальных показателей. 

Для достижения поставленной цели решены следующие задачи: 

— получение уравнений регрессии отношения сигнал/шум для мультипликатив- 
ного шума на значения мультифрактальных показателей и оценивание качества этих 
уравнений; 

— разработка метода оценивания параметра мультипликативного шума по мульти- 
фрактальному показателю; 

— исследование характеристик разработанного метода на тестовых изображениях. 

Оценку отношения сигнал/шум для мультипликативного шума на изображениях 
целесообразно проводить на основе зависимости статистических характеристик изо- 
бражения от масштаба данных. 


Изложение основного материала 


При анализе зависимости статистической информации от масштаба данных 
используется два подхода: структурные функции и сингулярные меры [8]. Под струк- 
турными функциями понимают статистические моменты абсолютных разностей, за- 
висящие от масштаба на некотором интервале его значений. Приведем определение 
структурной функции для одномерных данных, представляющих значения интенсив- 
ности 1-го столбца изображения /(х,у,), х=1,..., М ‚ где М- число пикселей в столбце 


изображения. 
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Разностью 7-го порядка функции [(х, у,) называется 
АЦг,х, у.) = Цх+г,у,)-ТОьу,), х=1...,М-г. 


Определим степенную зависимость с(7т) момента порядка т У ИЕ Зы от 


‚ й 5(т) 
}*[^) ‚ 120, () 


где (.) — означают усреднение по переменной х. 


масштаба г как 


(Аи, У) 


На основании (2) в работе [8] введен мультифрактальный показатель 
Н =5(}. 

Наибольшее возможное значение для Н! = 1 достигается для почти всюду диф- 
ференцируемых функций. Если Н! = 0, Г(х, у,) — стационарный процесс, график кото- 
рого заполняет пространство. 

На рис. 2 показан характер зависимости мультифрактального показателя Н! от 
вида помехи для столбца тестового изображения, содержащего один объект. Н\ в этих 
случаях равно 1,0133 и 0,99685 соответственно. 


1 1 
Вх) Их) 
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пространственная координата х, м пространственная координата Хх, м 
а) 6) 


Рисунок 2 — Значения интенсивности столбца тестового изображения (а) и они 
же для изображения, пораженного мультипликативной помехой (6) 


Для вычисления Н! в результате логарифмирования (2) получаем 
105 (| (",х,у,) й = 6(г)108 |) + С, 


где С — некоторая константа. Значение Н, = 6(Т) находим путем определения пара- 


") от к 


Другим подходом к мультифрактальному анализу данных являются сингулярные 
меры, представляющие собой локально усредненный модуль градиента изображения. 

В основе определения сингулярной меры лежит понятие нормализованной раз- 
ности 1-го порядка. Ее значения вычисляются в каждой точке столбца изображения 
по формуле 


метров линейной регрессионной зависимости 105(|АГ(г,х, у.) 
1 


Ах, у,) 
1 = Раны 
&( и) (АГЬх,у,) Хх 


где (1АГ(,х, у, = УХА») 


х= 
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Пусть масштаб г принимает значения 1, 2, 4, ..., №. Тогда сингулярная мера 
=(г,х, у, ) определяется как 


х+и-1 


а УЕ =, Мы 
Г хх 


Определим функцию К(т), выражающую степенную зависимость сингулярной 
меры от масштаба, по формуле 


—К (т) 
Г 
$. 
(ебьх,у,)г) «| — ‚ т>0. (3) 
№ 

В смысле сингулярных мер значения интенсивности изображения характери- 
зуются мультифрактальным показателем 


1 
С =К'(0>0. 
Наибольшее значение Ст, равное 1, достигается для случайно расположенных 
д -функций, вся энергия которых сконцентрирована в отдельных точках. Случай С, > 0 


(наименьшему значению) достигается для случайно расположенных функций Хевисайда. 
Таким образом мультифрактальный показатель С! характеризует уровень импульсного 
шума на изображении. 

Для вычисления С! в результате логарифмирования (3) имеем 


105 (Еб,х, у," =-А (т) [т Е", 


где С’ некоторая константа. Значение С! определяем как аппроксимацию производной 
функции К(т) в точке т =1. 

Приведенные определения структурной функции и сингулярной меры позволяют 
вычислять мультифрактальные показатели Н! и С! для одномерных данных. Для обра- 
ботки изображений, представленных матрицами дискретных отсчетов интенсивности, 
введены показатели Н\, С! по направлениям х, у: Н‚, Н‚, Ну», С», С» Су, соответственно. 
При определении Н,, С. разности 7-го порядка вычисляются и усредняются по строкам 
изображения. При вычислении Н,, С, обработка производится по столбцам изображения. 

В процессе определения Ну» С», разности т-го порядка вычисляются и усред- 
няются по двумерным окрестностям пикселей изображения. 

В работе [9] получены графики зависимостей значений мультифрактальных пока- 
зателей от отношения сигналлпум для мультипликативного диагонального шума (рис. 3) 
и вертикального (рис. 4) на тестовом изображении. 

У мультипликативного шума вертикального и мультипликативного горизонталь- 
ного размеры идентичны, т.к. изображение центрально-симметричное, а мультипликатив- 
ная помеха моделируется детерминированной функцией пространственной координаты. 

На основе этих зависимостей в данной работе осуществляется оценивание па- 
раметра мультипликативного шума. Получим уравнения регрессии, описывающие зави- 
симость отношения сигнал/шум 4 по мощности от значений показателей Н,, Н,, Ну», 
С,» С». В качестве регрессионной зависимости выбрана полиномиальная. Степень 
полинома и является параметром метода оценивания. 

Коэффициенты полинома регрессии оценивались методом наименьших квадратов, 
который используется для приближения функций, заданных числовым массивом [10]. 

Полученные регрессионные зависимости изображены для мультипликативного 
диагонального шума на рис. 5 и мультипликативного вертикального шума на рис. 6. 
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Рисунок 3 — График зависимости мультифрактальных показателей Н, (а), Н, (6), 
Нуу (в), С. (г), С, (д), Су» (е) от отношения сигнал/шум для мультипликативного 
диагонального шума 
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Рисунок 4 — График зависимости мультифрактальных показателей Н,.(а), С. (6), 
С (в) от отношения сигнал/шум для мультипликативного вертикального шума 
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Рисунок 5 — Регрессионная зависимость 4 по мощности 
для мультипликативного шума от Н, (а), Н, (6), Н.» (в), С, (г), С, (д) С,» (е) 
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Рисунок 6 — Регрессионная зависимость 4 по мощности 
для мультипликативного вертикального шума от Н, (а), С, (6), С», (в) 
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Качество уравнений регрессии отношения сигналЛпум для мультипликатив- 
ного диагонального шума на значения Н,, Н,, Ну», С», Суи для мультипликативного 
вертикального шума на значения Н,, Су, С‚, определялось с помощью коэффициента 
детерминации по формуле [10]: 


В =1-А 
ет 
5 
где И дисперсия остаточных данных, А, — дисперсия исходных данных. 


Значения коэффициентов детерминации для определения качества аппроксимации 
полиномом 8-й степени отношения сигнал/шум для мультипликативного диагонального 
от значений Н,, Н,, Н» С» С» Съ и вертикального шума от значений Н‚, С,, С» 
приведены в табл. 1. 


Таблица 1 — Коэффициенты детерминации при аппроксимации зависимости 
отношения сигнал/шум для мультипликативного диагонального шума от значений 
мультифрактальных показателей Н,, Н‚, Н,», С», С» Су 


Показатель Дисперсия Дисперсия исходных | Коэффициент 

остатка данных детерминации 

Н; 0,1057 1,3301е+003 1,2585е+004 

Н, 0,1057 1,3301е+003 1,2585е-+004 

Ну 0,2147 1,3301е+003 6,1961е-+003 

С» 2,1972 е-025 1,3301е+003 6,0538е+027 

С, 2,1972 е-025 1,3301е+003 6,0538е+027 

С 2,3979 е-023 1,3301е+003 5,547 1е-+025 
Таблица 2 — Коэффициенты детерминации при аппроксимации зависимости 
отношения сигналлпум для мультипликативного вертикального шума от значений 


мультифрактальных показателей Н,, С», Су 


Показатель Дисперсия Дисперсия исходных | Коэффициент 
остатка данных детерминации 
Нх 0,0065 1,3301 е+003 2,0475 е+005 
Е 9,5264 е-028 1,3301 е+003 1,3963 е+030 
Су 43,5515 1,3301 е+003 30,5419 


Из шести исследованных зависимостей выбрана та, для которой коэффициент 
детерминации максимален. Для каждого случая — это зависимость 4 от С,. Она 
аппроксимируется полиномом 8-й степени: 


а=а +асС, + О. + а, С’ +... + а, С”, (4) 


где а,, 4, ... ‚а, — коэффициенты аппроксимирующего полинома. 

С учетом этой формулы метод оценивания отношения сигнал/шум для мульти- 
пликативного диагонального и вертикального шума сформулируем следующим образом: 

— для исходного изображения оценить Сх; 

—с использованием формулы (4) вычислить отношение сигнал/шум для мульти- 
пликативного диагонального и вертикального шума. 
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Для оценивания Сх предложен следующий алгоритм. 

1.Задаем максимальный размер окрестности для усреднения градиентного 
ПОЛЯ Глах = 30 и задаем последовательность значений порядка момента 4 из интерва- 
ла (0,5; 1,5) с шагом 0,25. 

2. Для определения значения мультифрактального показателя Сх вычисляем 
градиентное поле изображения по формуле: 


[6.7 х, у) | = | (хУ+П- Г(х, у) 


3. Усредняем градиентное поле изображения по окрестности размера хи, 
изменяя 7 от | до Гиах. В результате получаем для каждого г изображение Б(х’ х, у) 
такого же размера, как исходное. 

4. Полученное изображение Р\(и, х, у) возводим в степень 4 поэлементно и 
вычисляем среднее абсолютное значение интенсивности т(г, 4). Получаем после- 


и 


довательность значений (г, 4), "=, ..., Гиах: 
5. Для каждого значения 4 выполняем линейную регрессию [и(и(х, 4)) на [и(х), 
Г=1,..., Гиах. Значения коэффициента при старшей степени уравнений регрессии для 


разных 4 формируют последовательность 5(4). 

6. Выполняем регрессию полиномом 4-й степени /[и(5(4)) на (4). 

7. С учетом значений коэффициентов полинома регрессии мультифракталь- 
ный показатель С, вычисляется по формуле: 


С, =4а4 + Заз + 2а2 +аь, 


где а1, 42, аз, ал — коэффициенты регрессионного полинома, соответствующие степе- 
ням 1,2, 3, 4 соответственно. 


Экспериментальные исследования и ВЫВОДЫ 


Экспериментальные исследования предложенного метода проведены на изобра- 
жении размером 256х256 пикселей, в центре которого — белый квадрат 64х64 пикселя 
на черном фоне. На это изображение накладывались следующие виды помех: мульти- 
пликативный горизонтальный, вертикальный, диагональный шум, для которого 
К(х, у) = Азш 2лажА ‚ где частота дискретизации Д = 1 / 256, в = 0,5, А определяется, 


чтобы отношение сигнал/шум 4 е {100; 50; 20; 10; 5;2; 1$}. 

Для полученных изображений вычисляем мультифрактальный показатель С,. 
Согласно предложенному алгоритму и оценивался параметр отношения сигнал/шум 
для мультипликативного шума по формуле (5). 

Оценка количества полученных результатов выполнялась с помощью относи- 
тельной и абсолютной погрешностей. Эти погрешности определяются следующим 
образом. 

Пусть а — точное значение, а* — приближенное значение некоторой величины. 
Абсолютной погрешностью приближенного значения а* называется величина 


А(а*) = |а —а ы . Относительной погрешностью значения а* (при а*0) называется 
А(а*) 


величина О(а*) = Е : 
а 


Результаты вычисления относительной и абсолютной погрешностей для оцени- 
вания отношения сигнал/шум для мультипликативного диагонального шума приведены 
в табл. 3 и мультипликативного вертикального шума в табл. 4. 
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Оценивание параметра мультипликативного шума на изображении... 


Таблица 3 — Результаты оценивания отношения сигнал/шум 


для мультипликативного диагонального шума предложенным методом 


ДЕ 


Истинное значение Оценка Абсолютная Относительная 

параметра шума параметра шума погрешность погрешность 
100 99,999 0,6253 е-012 0,0063 е-012 

50 49,999 0,1137 е-012 0,0023 е-012 

20 19,999 0,0568 е-012 0,0028 е-012 

10 9,9999 0,2274 е-012 0,0227 е-012 

5 4,9999 0,6253 е-012 0,1251 е-012 

2 1,9999 0,4547 е-012 0,2274 е-012 

1 0,9999 0,6821 е-012 0,6821 е-012 


Таблица 4 — Результаты оценивания отношения сигнал/шум 


для мультипликативного вертикального шума предложенным методом 


Истинное значение Оценка Абсолютная Относительная 

параметра шума параметра шума погрешность погрешность 
100 99,9999 0,1279 е-012 0,0013 е-012 

50 49,9999 0,0497 е-012 0,0010 е-012 

20 19,9999 0,0568 е-012 0,0028 е-012 

10 9,9999 0,0568 е-012 0,0057 е-012 

5 4,9999 0,0462 е-012 0,0092 е-012 

2 1,9999 0,0782 е-012 0,0391 е-012 

1 0,9999 0,1030 е-012 0,1030 е-012 


Сравнение полученных результатов оценивания параметра мультипликативного 
диагонального и вертикального шума предложенным методом (табл. 3, 4) прово- 
дилось с результатами оценивания этого же параметра методом работы [11] (табл. 5, 6). 
Для реализации метода работы [11] в зашумленном изображении считали стандартное 


отклонение $ и вычисляли отношение сигналлпум по формуле а= де 


т 
514? 
п = 0,0925 для диагонального шума ип = 0,1920 для вертикального шума. 


Таблица 5 — Результаты оценивания отношения сигнал/шум 


для мультипликативного диагонального шума методом работы [11] 


Истинное значение Оценка Абсолютная Относительная 

параметра шума параметра шума погрешность погрешность 
100 99,9993 0,7060 е-003 0,7060 е-005 

50 49,9996 0,3530 е-003 0,7060 е-005 

20 19,9999 0,1412 е-003 0,7060 е-005 

10 9,9999 0,0706 е-003 0,7060 е-005 

5 4,9999 0,0353 е-003 0,7060 е-005 

2 1,9999 0,0141 е-003 0,7060 е-005 

1 0,9999 0,0071 е-003 0,7060 е-005 


Анализируя полученные результаты, заметим, что для мультипликативного 
диагонального шума абсолютная и относительная погрешность оценивания предложен- 
ным методом меньше, чем методом работы [11] в 4 раза при значениях отношения 
сигналлпум 1 — 100 по мощности. Для мультипликативного вертикального шума абсо- 
лютная и относительная погрешность оценивания предложенным методом меньше, 
чем методом работы [11] в 3 раза при значениях отношения сигнал/шум 1 - 100 по 
мощности. 
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Таблица 6 — Результаты оценивания отношения сигнал/шум 
для мультипликативного вертикального шума методом работы [11] 


Истинное значение Оценка параметра Абсолютная Относительная 

параметра шума шума погрешность погрешность 

100 100,0512 0,0512 5,1205 е-004 

50 50,0256 0,0256 5,1205 е-004 

20 20,0102 0,0102 5,1205 е-004 

10 10,0051 0,0051 5,1205 е-004 

о 5,0026 0,0026 5,1205 е-004 

2 2,0010 0,0010 5,1205 е-004 

1 1,0005 0,0005 5,1205 е-004 
Таким образом в данной работе получены уравнения регрессии, которые позво- 
ляют оценить зависимость отношения сигналлпум тестового изображения, пораженного 


мультипликативным шумом, от значений мультифрактальных показателей. Разработан 
метод оценивания параметра мультипликативного шума по мультифрактальному 
показателю. Результаты исследований показали, что разработанный метод оценивания 
параметра мультипликативного шума целесообразно применять для выбора типа и 
параметра предварительной обработки. Направлением дальнейшей работы является 
разработка новых методов оценки параметров мультипликативного шума. 
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Ю.В. Емець 

Ощнювання параметра мультиплкативного шуму на зображенн! за допомогою 
мультифрактальних показник в 

Розроблено метод ощнювання параметра мультиплкативного шуму за мультфрактальним показником. 
Дослджен! характеристики розробленого методу на тестових зображеннях. 


Уи.Г. Ете5 

Мша@рНсабуе М№обе Рагатеег Еуашайоп оп Ппаре Бу шеап$ 07 Мш@йгасва| ш@дехе$ 

Тве тефо4 оЁ ша @рИсайуе по15е рагатеег еуалайоп аё а пи!Насва| шдех 15 деуеоре4. СБагасетзИс$ 
ое деуеоре4 тео4 аге ргоБе4 оп Фе {е5{ ипазез. 
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